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Abstract 

A new model to calculate the rate of nucleation is formulated. This 
model is based on the classical nucleation theory but considers also 
vapor depletion around the formed embryo. As the result the free 
energy has to be recalculated which brings a new expression for the 
nucleation rate. 

Introduction 

A typical example of the first order phase transitions is the condensation 
of a supersaturated vapor in a state of liquid droplets. It is wide spread in 
nature and became a model for theoretical investigations. Namely in this case 
the most advanced models of the phenomena were presented. The start of 
systematic investigations was given by the famous Wilson's chamber already 
at the end of the 19-th century pQ. Despite the long history there are still 
some unresolved problems in coincidence of the theoretical description with 
experimental results. Experimental technique is rather accurate, theoretical 
models are sophisticated but the measured rate of nucleation can differ many 
times from the theoretical predictions. 

Creation of the so-called "classical theory of nucleation" by Becker and 
Doering [2], Volmer and Weber [3], Zeldovich [I], Frenkel [5] gives the expres- 
sion for the stationary rate of nucleation which is rather evident but contains 
the free energy of the critical embryo as a parameter. The last value has 
to be derived in frames of some model for the critical droplet. The reliable 
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model is the capillary approximation which applies the standard conceptions 
of thermodynamic to the droplet. 

The classical theory of nucleation was mainly completed at 40-ties but 
experimental measurements did not confirm predictions of the theory - the 
experimental values of the rate of nucleation can many times differ from 
theoretically given values. The most radical reconsideration of the theoretical 
description was given by Lothe and Pound [8] who suggested that inside the 
droplet a special account of rotation degrees of freedom has to be made. This 
initiated a wide discussion, many arguments were presented but there is no 
strong arguments to state that all rotational degrees of freedom have to be 
accounted in a special way. The only evident modification is to make a special 
account of rotational degrees of freedom for the embryo as a whole object. It 
was done by H. Reiss, J. L. Katz, E.R. CohenJT] but this modification can not 
bring the theory in coincidence with experiment. Honestly speaking one has 
to declare that the theory of Lothe and Pound is also far from experimental 
results for some substances meanwhile it gives some promising results for the 
water condensation. After numerous strikes and discussions in 60-ies and 70- 
ies the problem of discrepancy between theoretical and experimental results 
has not yet disappeared. 

One has to stress that the difficulty is not in the form of the rate of 
nucleation which is 

Jtheor = N eq Zexp(-F c ) 

where F c is a free energy of the critical embryo, N eq is a normalizing factor 
of the equilibrium distribution, Z is so-called Zeldovich factor. The main 
technical efforts of the classical nucleation theory were intended to determine 
Z. 

But the problem is also not in determination of Z which is evidently pro- 
portional to kinetic coefficient W + (which is the inverse mean time between 
collisions of the critical embryo with molecules of vapor) and to the con- 
densation coefficient a c . The typical error in Z under the rough estimation 
results in some essential value, but not in order of this magnitude. So, the 
error which is the source of discrepancy between theoretical and experimen- 
tal results can be made in N eg or in F c . Since these values appear together 
in the value of the equilibrium distribution N eq exp(— F c ) / exp(l) which is a 
boundary condition at the left boundary of the near-critical region one can 
not separate them. Historically it is preferable to speak about the error in 
F c regarding N eq as some smooth value like the number of monomers in a 
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system. 



Formulation of the problem 

One can say that the theoretical models are not so sophisticated and their 
further development can bring the theory in coincidence with experiment. 
Creation and development of mult i- component theory show that the flow 
will go mainly through the saddle point and the free energy F s at the saddle 
point can be only lower than the free energy of the critical cluster calculated 
in the one-component theory. The reasons for the multi-component con- 
sideration can be very various, for example, the existence of some compact 
complexes of molecules inside the embryos. But it is clear that any more 
sophisticated consideration can lead only to the decrease of the free energy 
of the critical embryo and, thus, to the increase of the nucleation rate. But 
the the experimental rate of nucleation J exp can be less than the theoretical 
rate J them- for some substances and the reconsideration of the theory can only 
diminish the last value. So, in the situation 

Jtheor 

the new more sophisticated models of the embryo can only enlarge the dis- 
crepancy. 

Certainly, one can attribute the discrepancy between theoretical and ex- 
perimental values to the coefficient of condensation a c . But it is hard to 
believe that the coefficient of condensation can be like 10~ 5 ^~ 7 which is nec- 
essary to ensure the coincidence. This value does not correspond to the 
values 10*0 -. — 2 of the condensation coefficient observed in experiments for 
the condensation in the bulk liquid. Really, the small droplet and the bulk 
liquid are different objects but one and the same molecular content of these 
objects makes it difficult to believe in such difference. 

Another possibility is to attribute this discrepancy to the microscopic 
corrections for the free energy. These corrections are very complex functions 
and are quite different for different substances. But here appeared another 
problem. The characteristic weight of these corrections for big embryos has 
to be relatively small. So, there has to be a sharp functional dependence of 
the ratio of the rates without and with essential corrections on the size of the 
critical embryo. But the real situation shows a smooth value for the ratio 
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of theoretical and experimental rates of nucleation. There is no evidence why 
in decomposition on inverse radius of embryo for microscopic corrections the 
third term has to attain an enormously big value. 

As the result we see that there has to be some another mechanism which 
is responsible for the decrease of the nucleation rate in comparison with the 
classical consideration. One of the possible answers is that the free energy 
which stands in the main term for the theoretical rate 

Jtheor ~ exp(-F c ) 

is not a true equilibrium free energy. The minimal work of the embryo 
formation can be attained only at the infinitely slow process of the embryo 
formation and here the situation is not so slow as it is necessary to apply the 
classical expression for the free energy 

F{v) = av 2/3 - bv 

of the embryo with v molecules where a is the normalized surface tension 
and b is the difference of the bulk chemical potentials. 

Here a careful reader has to argue because in the classical theory of nucle- 
ation there exists a solid justification to apply the thermodynamic expression 
for the free energy. We shall repeat here this justification. 

The embryos are so small and the liquid is so dense in comparison with 
vapor that the time of the diffusive and thermal relaxation is very small in 
comparison with the mean time between sequential collisions of the embryo 
with molecules of substance in vapor. The time of relaxation can be estimated 
as 

t rel = R 2 JAD 1 

where R d is the radius of the embryo and D h is the coefficient of diffusion in 
liquid. The time between sequential collisions is W+ 1 where W + is the direct 
flow of the molecules on the embryo. In the free molecular regime it can be 
estimated as 

W + pa -v t nS 

where v t is the mean thermal velocity of molecules in vapor, n is the density 
of the particles number, S is the square of the embryo. One can check the 
strong inequality 

trel < W" 1 
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and see that it really takes place. 

The point of view adopted in this theory states that there is no strict 
connection between this justification and the possibility to write the free 
energy only for the embryo. The embryo at such big times as the time of the 
while evolution can perturb the environment meanwhile it is a quasi-isolated 
object at the relatively small times. All adiabatic approaches deal with such 
situations. But nobody states that it is necessary to forbid the change with 
environment and to use only the static characteristics. 

So, here we shall consider the environment of the embryo and require that 
the embryo has to be in the stable state being embedded in this environment. 

Equilibrium environment 

There exists one argument which is rather unpleasant for the formulated non- 
stationary theories. Really, after the precritical embryo is dissolved the vapor 
environment is still perturbed. Now there is an excess of the molecules in 
vapor near the place where the embryo was previously formed. So, the story 
of this fluctuation continues. In some cases it will be completely dissolved, 
in some cases it will give a new life for a new embryo. So, we see that we 
come to a level of fluctuations in vapor. We need to consider not a single 
embryo but a fluctuation of molecules density where there is an embryo and 
this embryo is in equilibrium conditions. 

The main supposition used here will be a requirement to build the fluctu- 
ation of a minimal size corresponding to a possibility to have an equilibrium 
embryo. Later this supposition will be explained and justified. 

Construction of the theory 

Considering the fluctuation containing the liquid phase and the vapor phase 
we need to calculate the minimal work to form such a fluctuation. This work 
has to correspond to the equilibrium process and, thus, the fluctuation has 
to be an equilibrium one. Then, we turn to the question what fluctuation 
will be equilibrium one. 

Consider the system of N molecules and an embryo of v molecules inside 
it. At this very moment the number N is unknown but later we shall express 
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it through v. The substance balance gives 

N = v + q 

where q is the number of molecules in vapor. The density of the number of 
molecules in vapor phase is 

n = q/V 

where V is the volume of the fluctuation. If the initial density n of the vapor 
molecules is known then 

V = N/n 

Then in the expression F = av 2 ^ — bv one has to recalculate the difference 
in chemical potentials b on the base of n instead of n . Namely, we have to 
substitute n instead of n in the expression 

b = ln(n /noo) 

where is the molecular number density in the saturated vapor. Earlier 
this expression was 

6 = ln(C + l) C = (n - noo)/noo 

now it is 

b = ln((n - nooVnoo + 1) = ln(((7V - u)n /N - TiooV^oo + 1) 
Then 

F = av 2/3 - z/ln(C + 1) - u\n(l - ^) 

Since the fluctuation has to be the equilibrium one to ensure the minimal 
work of formation we seek the local minimum of the free energy F and the 
argument of this minimum has to be the true number v e of the molecules 
inside the embryo 

u e = arg(min u F(i / , N)) 

So, for given u e we get the value of N which provides the minimum of F 
namely at v e . It gives the dependence iV = N(u e ). The value v e is the true 
value of v to be used if we want to calculate the free energy of the embryo in 
the stable fluctuation. So, we get then the value F e as the value of the local 
minimum of the free energy, i.e. 

F e = F(v e ,N(v)) 
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This value will be the true value of the free energy of the embryo containing 
v molecules in frames of this model which will be called "the model with 
the minimal environment". We shall keep the name F eq for the free energy 
calculated in this model 

Expression for dFjdv looks like 

dF/du = \av~^ - ln(C + 1) - In(l - -£) + (1 

and the equation of extremum 

^- 1 /3 = ln(C+1)+ln(1 _^ ) 



N> (1-%)N 

is very simple but can not be solved analytically. So, it is not clear whether 
the free energy F(u, N) really have necessary minimum (the minimum at the 
origin does not match). So, it is necessary to see the properties of the free 
energy in this model. 

Instead of analyzing the dependence N{y e ) we shall consider dependence 
v e (N) and we see immediately that when N goes to infinity then v e also goes 
to infinity and F e will attain a big negative value. This shows that if for some 
v we can ensure the minimum, i.e. v = v e then one can provide minimum 
for all v greater than the mentioned one. The problem is for small v one can 
not get N providing v e . One can see it from considerations below. 

We shall only diminish the action of correction term if we linearize the 
logarithm. It is so because v/N is positive and In is a convex function. Then 

F = az/ 2 / 3 -z/ln(C + l) + ^ 
For this function one can easily calculate the second derivative 

d 2 F/du 2 = -(2a/9K 4 / 3 + ^ 

Then for the transition to the absence of the threshold we have 

v = (iVa/9) 3/4 

We see that the critical (minimal) value for iV ensuring a threshold is 

N c = 9z/ 4/3 /a 
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which is a big value and the argument of decomposition v/N is small at v >> 1 
which a natural requirement for thermodynamics. Then the decomposition 
is quite possible. 

For the critical value of N c we have 



9^ 4 /3 

or 

F — a(l + l/9)z/ 2/3 -is\n(( + 1) 

Really, we need only the barrier value of the free energy which corresponds 
to N = N c . 

So, we come to the situation with a redefined surface tension 

a ->• a(l + 1/9) 

The surface tension grows 1 + 1/9 times. 

This conclusion is very important. On one hand it allows to formulate a 
new variant of the Gibb's theorem about the proportionality of the free energy 
of the critical cluster to the part of the surface energy. Now all constructions 
will be the same but only the quantity of one third will be redefined. On 
the other hand it explains numerous conclusions in experimental papers that 
everything would be good with theoretical predictions if the surface tensions 
is redefined. Now the reason is obvious. 



The justification of the model 

Now it is clear why it was necessary to require the fastest creation of the 
critical embryo or the fastest transition over the activation barrier. We see 
that the precritical cluster is very unstable formation (there is even no stable 
environment). When the molecules will be spread the volume with the linear 
size like 

R w = (wv) 1/3 

there will be no trace of the embryo. Here v v is the characteristic volume 
for one molecule in the saturated vapor. The characteristic time can be 
estimated as 

t w = R 2 J(AD) 

where D is a coefficient of diffusion in vapor. 
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Here we take the longest estimate corresponding to the diffusion relax- 
ation while the disappearance of the embryo can occur faster. 

We have seen that to make the stable configuration it is necessary to 
spread the volume of the influence up to the linear size 

R q = (N c v v y/ 3 

One can see that 

Rq/Rw ~ v Xl% > 1 
The characteristic time here can be estimated as 

t q = R 2 q /(W) 

Since 

tq/t w ~ V 2/9 > 1 

it is necessary to gather the embryo as soon as possible and, thus, to gather 
the molecules from the minimal volume ensuring the stability of the fluctu- 
ation. Namely this was done in this model. 

The rate of nucleation 

Since Z does not strongly depend on supersaturation one can take it in 
expression for the rate of nucleation being calculated in frames of the old 
classical theory. This helps to ignore the absence of explicit barrier at the 
situation with minimal environment. The value of the normalizing factor 
N eq also can be taken from the old variant of the nucleation theory. One has 
to stress that here one can use as the reference theory all already existing 
approaches including Lothe and Pound approach [6], Reiss, Kats, Cohen 
approach [7J, etc. All already formulated approaches have nothing in common 
with depletion and the finite time of formation of the critical embryo and 
deals with another effects. So, they can be really combined with this approach 
without any difficulties in justification. Technical difficulties are also absent. 

One has to stress that here it is possible to avoid to linearize the logarithm 
of the new (correction) term in the expression for the free energy. Expressions 
will be slightly more complex but this approach will give the better accuracy. 
This improvement is important because the small relative shift in the height 
of barrier will radically change the value of the nucleation rate. Since the 
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formulas without this simplification are not explicit ones it is preferably to 
use them directly in numerical calculations. The scheme of approach remains 
the same. 

Here the main attention is payed to the period of formation of the em- 
bryos. It can seem that this period is not too important because in the 
process of the global condensation one can see the universal features [9], |10j . 

The final asymptotic of over-condensation (coalescense) is universal [TS] 
also. But the recent investigations of the over-condensation [8], [12] show that 
the initial data are very important in the final stages of the process and can 
determine the final behavior even qualitatively. 
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